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Abstract 

Multistate generalizations of Landau-Zener model are studied by summing entire series of per- 
turbation theory. A new technique for analysis of the series is developed. Analytical expressions 
for probabilities of survival at the diabatic potential curves with extreme slope are proved. Degen- 
erate situations are considered when there are several potential curves with extreme slope. New 
expressions for some state-to-state transition probabilities are derived in degenerate cases. 



PACS numbers: 75.10.Jm, 03.65.-w 



I. INTRODUCTION 

The famous Landau-Zener two-state model, introduced and solved in 1932 by Landau ^ , 
Zener {2], Majorana and Stiickelberg {4] finds many applications in atomic physics and 
beyond. This is due to its virtue of describing generic case of non-adiabatic transitions in 
quantum mechanics. The main feature of the exactly solvable quantum model is the linear 
dependence of the matrix Hamiltonian on time. The model allows the derivation of exact 
expression for the state-to-state transition probability. 

The natural generalization of the two-state model is the model with arbitrary (but still 
finite) number of states, A^. The linear dependence of matrix Hamiltonian on time is retained 
H(t) = A-|-Bt, where A and B are time- independent N x N matrices. Hereafter we show by 
bold type the operators and vectors in A^-dimensional linear space. The lower case characters 
are used to denote vectors in this space while the capital characters denote matrix operators. 

Without loss of generality one might assume that the basis is chosen in such a way that 
the Hermitian matrix B is diagonal, Bjk = PjSjk, where the parameters P have the meaning 
of slopes of linear diabatic potential curves. The so chosen basis is known as the diabatic 
basis. The non-diagonal elements of matrix A have the meaning of couplings between the 
diabatic states, Ajk = Vjk- The diagonal matrix elements of A play a different role. It is 
convenient to introduce for them a special notation, Sj = Ajj (these notations are the same 



as in our preceding studies ll|, ll3|, llJ]). The diagonal matrix elements of the Hamiltonian 



H(t) are referred to as diabatic potential curves. In the case of the multistate Landau-Zener 



f( 

The problem is to solve the non-stationary Schrodinger equation 



model, they are linear in time, E'^^^{t) = (3jt + Ej. 



= m c , (1.1) 

and to find S-matrix. Generally speaking the full exact solution of (11. ip is not available. The 
known exact solutions 

mm 

refer to special choices of the model parameters Sj, Vjk, 
such that the quantum interference oscillations do not appear in the transition probabilities. 
Furthermore, even in the case of the most general form of matrix Hamiltonian one can 
exactly find two elements of S-matrix which correspond to survival on the diabatic potential 
curves with extremal (maximum or minimum) slopes. The simple formula for such elements 

n 

was originally guessed by S. Brundobler and V. Elzer [3] based on numerical calculations. 
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The proof of Brundobler-Elzer (BE) formula was carried out recently by several different 
ways. A. V. Shytov obtained this formula via treatment within the contour integration 
approach 10| . M. V. Volkov and V N Ostrovsky carried out the proof using non-stationary 
perturbation theory [11]. However there are some oversights in this proof, as B. E. Dobrescu 
and N. A. Sinitsyn indicated in the comment to this paper p^]. The comment contains a 
new proof of BE formula partly based on developments by Volkov and Ostrovsky; at the 
crucial step it essentially uses results for the bow-tie model jo] exactly solved by Ostrovsky 
and Nakamura. 

The objective of the present study is to provide a proof of the BE formula which is devoid 
of deficiency of the previously suggested proof being fully based on analysis of non-stationary 
perturbation theory and summation of an entire perturbative expansion. Compared to the 
case of a Hamiltonian Hbound(^) with all the matrix element bounded [(-f^bound)jfc (t) < a for 
all times t] the case of the multistate Landau-Zener Hamiltonian provides important specifics. 
The emerging integrals typically contain highly oscillating exponential factors that ensure 
integral convergence. For some choice of parameters in the integrand the oscillations vanish 
which means that the integral is a singular function of parameters. These singularities are 
to be treated in the analysis with proper care; albeit namely the presence of singularities 
allows a closed-form evaluation for each term of the entire series with subsequent analytical 
summation. 

In the main Section UTTl we develop a new approach to treat the singularities. The prelim- 
inary Section [III introduces notations and contains a general description of the perturbative 
series. In distinction to the scheme suggested by Dobrescu and Sinitsyn [l2|], our proof (Sec- 
tion [TTll) does not use results of any exactly solvable model. We beheve that such a complete 
treatment of the perturbative expansion with analytical summation of series is of general 
interest . 

Another goal of our study is to consider some degenerate cases (Section llVj) . Here more 
state-to-state probabilities can be evaluated, up to the fully degenerate multistate model 
where an entire matrix of state-to-state transition probabilities is found (Section lIVDI) . 
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II. NON-STATIONARY PERTURBATION THEORY 



The well-known formula for transition probability for two-state linear model was derived 
by Zener by reducing the Schrodinger equation to an equation for a hyperbolic cylinder 

n fi 

function [2|. Majorana [3| used contour integration method in a complex plane to solve 
the same problem. Much later Kayanuma suggested an alternative approach [isl, [l^ where 
the non-stationary perturbation theory is used. As discussed in the Introduction, we in the 
present paper provide a generalization of this method to the multistate case. 

The non-stationary Schrodinger equation (11. ip might be written as the set of N coupled 
first-order differential equations: 



.dcj 



EjCj + PjtCj + ^ VjkCk 



j, k = l,2,...N . 



(2.1) 



After phase transformation which eliminates the diagonal elements on the right-hand side 
of equations (12.11] it takes the form 



.daj 
'~dt 



^ Vjk exp 



i { {e, - ek)t + -{(^j - (3k)t' 



J, fc = l,2,...iV. (2.2) 



The integral form of this equation 

aj{t) = aj(— oo) — i dt^ Vj^ exp 

k^i 



i{ej - ek)to + -{(3j - (3k)t, 



is convenient for an iterative solution. The successive approximations, a^"'\tn), are found 



(2.3) 



by iterations: 



(0), 



dtn ^ Vjk exp 

ky^j 



i{ej - ek)tn + ^^(/?i - I3k)tl 



We use label 1 for the initially populated state, so that initial populations aj{—oo) are 



aj(— oo) = 5ji . 



(2.5) 



Then transition probability to j-th state is 



lim a^"^(-|-oo) 



(2.6) 



In the next formula we introduce a vector-function of time f(t) = {fi{t), f2{t) . . . fN{t)}, 
which is a vector in A^-dimensional linear space. The operator T is x matrix; it 



transforms the vector- function f(t) into another vector-function with components: 



N 



Tf 



.(Wi) = {-i)^yjk / dtnexp 



/jk 



tn + l 



i{ej - ek)tn + - {13 j - Pk)tn 



fk{Q . (2.7) 



With respect to the time variable the operator T is an integral operator. Our equations 
(12. 4p can be written as 



(2.8) 



where dependence on time is implicit. The zero iteration a^^^ is defined by the initial 
conditions (12. 5p : a^°^ = 5ji. 

We further introduce the vector in A^-dimensional linear space by a formula describing 
its components dj: 



df\t) = -iVji I dtiexp 



(2.9) 



The j = 1 component d^i^ is assumed to be zero by definition. Similarly, the vector d^"''' 
(m > 2) is given as 



N 



N 



N 



N 



km — 2 -/-km — 1 



fci#i 



ft rtm 
X / dt^ / dt^^x ... dti 



— oo 



oo 



X exp 



X exp 



— oo 
m—l 



K^j - ^km-i)tm + i Y^^''' ~ + ii^k, - ei)t 

i=2 
. m—l 



1=2 



(2.10) 



If the couplings are small, then the order of magnitude estimates are T ~ d^"^) ~ V^. 
Note the important relations between operator T and vectors d^"^^: 



, ^, . . . , 



Xd("^) = d("+^) , m = 1, 2 
ta(o) = d(i) . 

Using these relations and equation (12. 8p we express the n-th iteration to a as 



f2.111 



m=l 



(m) 



(2.12) 
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Formula (12.120 is the basis for all subsequent analysis. In order to find some transition 
amplitude one should evaluate the corrections fl2.10p to all orders m in the limit t +00, 
then sum up all corrections using equation fl2.12p with n — > +00. The sought for probability 
is given by formula fl2.6p . 



III. PROOF OF THE BRUNDOBLER-ELZER FORMULA 
A. Preliminary transformations: change of variables 

Consider the case when the initially populated non-degenerated diabatic potential curve 
has ES, i.e. its slope is the largest = max^ jSj) or the smallest {Pi = min^ (3j) of all slopes. 
Here we set out to find the survival probability on such a potential curve. The general vector 
formula fl2.12p for the first component reads 



(m) 



(3.1) 



m=l 



m=l 



The arbitrary term in the sum is given by (12.100 and (12. 9p . In the limit t ^ 00 and after 
reducing the brackets we obtain: 

N N N N 



^m — 2 7^^m— 1 
ft2 



fcl#l 



r>00 

X / rftm. / titm—l 

-00 J —00 J —00 

m— 1 



X exp 



X exp 



dtm-l ■ I dti 

i=l 
. m—1 



i=l 



(3.2) 



Let us now introduce new integration variables {xi, . . . , Xm} such that: 



X- 
X 



m 1 



Xm. e -00, 00 



Xj G (0, 00), j = 1, 2 . . . , m — 1 . 



(3.3) 



The important advantage of this transformation is that the ranges of variation of the new 
variables are simple and unambiguous, cf . discussion in Refs. jl2 . fisl ] . The Jacobian of the 
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transformation is equal to (—1)™ ^, the inverse transformation is given by 



m— 1 



tj Xm ^ ^ ^A; ) j 1)2,..., TJl 1 , 

k=j 



(3.4) 



In order to express the integrand in ( 13. 2p in new variables the following formulas are useful: 



m—l 



A;=l 



»7l— 1 



''»7l— 1 



k=l 



^k=l 



m—l 



fc=j+l 



+2 

"m—l m 



z = 1, . . . , m — 2 ; 



(3.5) 



In new variables the integral is cast as 

N N 



N 



N 



- (-i)"" ^ Vlkrr^-i V'fc^-ifc^-a • • • Hafe VkiklVki 

r*©© /'oo /'oo r m—l m—l 

X / (ixr,!, / dxm-i ■ ■ ■ / c^a^i exp 



m 

oo ^0 







^^1 ^ ^ -^n ^ ^ ^ ^kn-^n 
n=l 



n=l 



X exp 



X exp 



m—l 



m—l 



A 2xm - 



n=l 



, n=l 



m-2 



m—l 



^ ^ j 2xm( x„) -|- j -|- 2 ^ ^ 



ra=l 



i=n+l 



X exp [{-2x^Xm^i + X^_i) 



(3.6) 
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The integration over dxm in infinite limits gives a (5-function. After reducing brackets in tlie 
exponents one obtains 

N N 





N 




E 


km-lj^l 


km — 2 -/-km — \ 


POO 


POO 




/ dxi exp 


Jo J 






m— 1 



k2j^k3 



^knjXr 



ki^k2 



X exp 



x2n6 



. m—l 

I 

2 



n=l 
m-2 



m—l 



n=l 



n=l 



j=n+l 



m—l 



E(^l ~ f^kJXr^ 



n=l 



(3.7) 



Tlie subsequent analysis of the multiple integral in (13.71) essentially depends on how much of 
the indices fc„ are equal unity. At first we will consider the case when all indices are different 
from unity. Subsequently the integral with an arbitrary set of indices will be evaluated. Note 
that in this subsection we do not use the condition that the slope (3i is extremal. However 
in the next subsection this assumption becomes essential. 



B. The case with kn ^ I for all n 



We carry out a new change of integration variables in such a way that the argument of 
the ^-function in (13.71) depends on a single new variable: 



Vi = - Pkjxn , i = 1, 2, . . . , m - 1 . 



(3.8) 



n=l 



The integration limits in the new variables has a simple form due to the fact that Pi has 
extreme value compared with all another slopes. For the sake of definiteness we assume that 
Pi = maxj /3j, then 



Vm-i G (0,oo) , 
Vi e (0,2/j+i), 



, m — 2 



(3.9) 



The modulus of the Jacobian for this transformation is 

m—l 



(3.10) 



Let us denote the multiple integral in (13.71) as /. Then in new variables we have 

I = 2ir\J\ dym-i6{ym-i) j dym-2--- rfyi /(yi, ?/2, • • • , ^m-l) , (3.11) 

Jo Jo Jo 

where f{yi, y2, ■ ■ ■ , ym-2, Vm-i) is a regular (smooth) function of all its arguments. One can 
see that the integration over dym-i with 6{ym-i) in the integrand implies that ym-i 0. 
This contracts the integration range over all other variables to zero. Thus, the entire integral 
/ is zero. 



C. The case with arbitrary set of indices 



Let us assume that {p — 1) of indices in (13.71) are equal to one, where p < m. Taking 
into account the obvious restrictions {ki ^ 1 and k^-i 7^ 1 and ki^i 7^ ki), one obtains a 
limitation for p: p < |m for even m and p < |(m — 1) for odd m. In order to evaluate / in 
this case we need new notations. Let us introduce a string of integers S = {si, S2, ■ ■ ■ , Sp_i} 
that includes all the labels Sj such that kg^ = 1. It is an ordered set, so that Sj+i > Sj. The 
complementary string C = {ci, C2, . . . , Cm-p} includes all labels cj such that kc^ 7^ 1 and also 
is ordered: Cj+i > Cj. The multiple integral in (13.71) is 



dXr 



exp 



neC 



X exp 



n&C 



nec 

n^m — l 



.n&C 



X / dxs^ / dxs2 ■ ■ ■ I dxsp_^ exp 
Jo Jo Jo 



n^C j>n 



(3.12) 



The integration variables belonging to S string enter exponent linearly (while other variables 
provide quadratic terms as well). This allows us to carry out integration in semiinfinite 
interval using the formula: 



00 1 

jfcx — 



e'^'^dk = iV- + -kS^x) 
X 



(3.13) 



9 



Here indicates integration in the principal value sense. After this (13.121) reduces to 



dx. exp 



X exp 



nec 



27t6 



n<j 



\ 



nec 



+ iV- 



1 



n<j 

- E iPl-PkJXr. 



Now a change of variables (13.81) is conveniently modified to 



(3.14) 



Vi = "^iPi - I3k,jxc„, i = l,2,...,m-p , 

n=l 

Vm-p e (0, oo) , 

e (0, Vi+i) , i = 1, 2, . . . , m - p - 1 . 



(3.15) 



The Jacobian modulus is 



(3.16) 



Each of delta-functions in formula (I3.14p depend only on single new variable yi, so that this 
formula is cast as 



oo rVm-p 

p ^\ym—p) I dyui—p—i 
^0 



/ = 2n\J\ / dym-p S{y 
Jo 

P-l r 



y2 



dyif{yi,y2,---,ym-p) 



TTSi-ys,-j) + iV- 



JSj-3 



(3.17) 



where f{yi,y2,---, ym-p) is a regular function of all its arguments. As in previous subsection, 
the integration over dym-p with 5-function contracts to one point, namely zero, the range 
of integration over all other variables; thus it could be said that the contribution from the 
P-terms is zero because of identity 



y I 

V-f{x) dx^O 
^ 



(3.18) 
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for y — ^ and f{x) non-singular at x = 0. Therefore the entire integral is different from zero 
only if integrand is singular function of all its variable. It could be only if the number of 
integrals in (13.1 7p equals the number of (5-functions in integrand. This reasoning give us the 
condition m — p = p — 1 + 1, i.e. m = 2p. This means that only even terms in the expansion 
fl3.ip give non-zero contributions. The string S consists of (|m — l) numbers. Taking into 
account the inequalities fcj+i ^ ki, ki ^ 1, fcm^i 7^ 1 we obtain the necessary condition for 
indices in (13. 7p : 

k2j = l for j = l, 2, ...,^m-l. (3.19) 
In other words the following indices have the value 1: 

^2) ^4) kQ, . . . km-4, km-2 ■ (3.20) 



D. Summation of non-zero contributions 

For an arbitrary term in (13. ip we obtain 

N N N 

P i>co r-Dp i>y2 

^ n — ^ — I / / ^yp-^ " 

J=l ~ Pfc2j-il Jo Jo Jo 

p 

xfiy,,y2,...,yp)l[S{y,) , p= 1, 2, .... (3.21) 

1=1 

The product of 5-functions in the last expression makes the integrand to be symmetrical 
function with respect to arbitrary permutations of the integration variables {yi, y2, ■ ■ ■ , yp}. 
Besides this, the integrand is an even function of any of its argument that allows us to 
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extend the limits of integration: 

N N N 



2p\ 



P POO roo POO y 

j—l IPl Pk2j-i\ J —oo J ~oo J—oo 

fe„-i^l fci^l .7=1 ^^=2,-1 1 



Here we used the property /(O, 0, . . . , 0) = 1. 

For the survival amplitude in the limit n — >^ oo we have the exact expression: 

N \ P / N 



-^E^ Ei^NexpUi:^), ,3.3, 



p=l ' \k^l "'^ /-'-I y y 

Finally, for the survival probability we obtain BE formula: 

Pn^M^)f^...[-2.t^^. (3.24) 

IV. EXTENSION OF THE APPROACH TO DIFFERENT DEGENERATE CASES 

In this section we assume the presence of a special property of a Hamiltonian compared 
to general treatment of previous section. Namely, we presume degeneracy of the potential 
curves. As above we consider the situation when the initially populated state 1 has the 
largest = maxj f3j) or the smallest (/^i = miuj f]j) of all slopes, except slopes for the 
states 1, 2 . . . , n (j 7^ 1, 2, . . . n). In other words, we presume degeneracy of extreme slopes, 
A = = ■ ■ ■ = Pn, or, in yet other words, there are n parallel curves with extreme slope. 
It is natural also to presume that the parallel curves are not coupled, i.e. Vij = for 
= 1,2, . . . 

The particular case when two bands of parallel potential curves received some attention 
in the literature 
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Subsequently we consider yet more special situation that the extreme slope curves are 
not only parallel, but fully degenerate, i.e. ei = £2 = ■ ■ ■ = £n- 
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A. The case of parallel diabatic potential curves with extremal slope 



In this subsection we consider the case of n diabatic potential curves with the same 
extreme slope (3i = [3 {i = 1,2, ... ,n) and (3 = maxfc>„{/?fc} or /5 = minfc>„{/?fc}. We also 
assume that 7^ Ej and Vij = 0, where (z 7^ j) and = 1,2, ... ,n). Such a model for 
n = 2 was considered in our previous work [14], where transition probablity Pu for (£2 > £1) 
was considered; now we concentrate on the survival probability. We will prove the formula 
for survival probability on the diabatic potential curve with extremal slope for n = 2. The 
proof for arbitrary n might be carry out similarly. 

The survival amplitude is again given by general formulas (13. ip and fl3.7p but the 
subsequent analysis is a little more complicated. The string of integers S is introduced as 
in the previous section. Besides this, we introduce a string of integers 71 = {ri, r2, . . . , Vg}, 
which includes all labels such that kr^ =2. It is also an ordered set: r^+i > r,. The 
complementary string C = {ci, C2, . . . , Cm-p-g} includes all labels Cj such that 7^ 1, 2 and 
also is ordered: q+i > Cj. The dimensions of these strings have to satisfy the conditions: 



p + g < —m 



for even m, 



P + g < 2^"^" ^) 



for odd m , 



(4.1) 



otherwise one or more of the couplings in (13. 7p is zero. The multiple integral in (13.70 is in 
this case after integration, given through the formula (13.130 : 



dXr 



dx 



C2 



dXr 



exp 



n&C 



X exp 



i65 



n<j 



nec 



jec nec 



2ti5 



n<j 
n&C 



n<j 



n<j 



(4.2) 



We realize the change of variables by formula (I3.15p with the same Jacobian modulus 
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fl3.16p . but now the total amount of variables is {m — p — g). Note that every delta-function 
after such transformation depends on only one variable. In new variables the multiple 
integral is given by the expression: 

I = 27r|J| / dy^^p^g 5{ym-p-g) / dym-p-g-l- ■ ■ / I/2, ■ ■ ■ , 

Jo Jo Jo 

i=i 

(4.3) 



xn 



-^H-Vsj-j-a,) + iV-_ 

Tr6{-yr^-j-fs^ +£1- £2) + iV- 



~yrj-j-l3j + £1 — £2 

Here aj is the number of the elements of the string TZ which are less than Sj, (3j is the 
number of the elements of string S which are less than rj and f{yi,y2,---,ym-p-g) is a 
regular function of all its arguments. Note that all 5-functions in the integral depend on 
different variables. 

The integration over dym-p-g with 5-function contracts to one point, namely zero, the 
range of integration over all other variables. Thus it could be said that the contribution 
from the P-terms is zero because of identity fIS.lSp . Furthermore the contribution from 
5-functions in the second product in (14.31) is zero. The multiple integral is different from 
zero only if integrand is singular function of every integration variable. This only happens 
if the number of integrals in (14. 3 p equals the number of 5-functions in the integrand. This 
reasoning give us the condition m — p — g = p — 1 + 1, i.e. m = 2p + g. Note that if 7^ this 
condition contradicts (14. ip . Thus, this implies that g = 0. Thereby we come to the same 
result: m = 2p as in previous section. Besides this, we obtain the complementary condition 
kj 7^ 2 for j = 1,2,. . .,m- 1. 

The same calculations as in non-degenerate case give us the survival probability 

PnH.Mr^e.p(-2.f Jl^). (4.4) 

For more general case of n-fold degeneracy {n < N) of extreme slope potential curves we 
similarly obtain 

P,-, = |ai(oo)|' = expf-27r V -^^^| j = l,2...,n. (4.5) 

In case when a band of parallel potential curves is crossed by a single curve {n = N — 1) 
formula (14. 5p reproduces an early result by Demkov and Osherov 

14 



B. The case of merged diabatic potential curves with extremal slope 



Consider the case when we have n diabatic potential curves with the same slope (3i = (3 
{i = l,2,...,n) and (3 = maxfc>„{/5fc} or (3 = miukyniPk} ■ As distinct from previous 
subsection we assume that Ei = e {i = 1,2, .. . ,n). This means that the potential curves 
1,2, ... ,n are merged. At first we will obtain expressions for survival probabilities for n = 2 
and then will generalize them for arbitrary n. 

In the case of two merged diabatic curves with extremal slope we assume the following 
conditions for couplings 

V2^ = C^Vu (4.6) 

with some z-independent constant C2. Acting further as in non-degenerate case we obtain 
restrictions for the coefficients: 

hj = l,2 for 3 = 1,2,..., \m-l. (4.7) 

For an arbitrary term in (13.11) we have after integrating: 

fc2p-l^l,2 fc2p-2 = l k2 = l fcl7^1,2 



P 

X 

J 

j2p-l 



p I 

niA-/3fc2.-j ' 



df-' = . (4.8) 

Due to the property (14.71) . summations over two terms, emerge here. Now we use 
condition (14. 6 P to get: 

2 

^k2j-lk2jVk2jk2j + l = Vfc2,„ilVlfe2^ + , + Vfc2^_i2^2fc2j + 1 = + ^2) Vk2j-llVlk2j + l ■ (4-9) 

k2j=l 

Then formula (14. Sp is rewritten as 
1 + cU ^ 

2 \fc5^1,2 



;i + 4)nv,kVki ] ]_ 

..,1,2 1/^^-/5^1 J 



df^' = . (4.10) 

15 



Obviously, are terms in the expansion of an exponent, 



[1 + 4)7rVM i 

lA-AI 



(4.11) 



However, the first term in formula (13. ip is 1, that is different from the first term in the 
expansion of expression (14. lip . This is easily taken into account. For survival probability 
we thus obtain 

2 



Pn 



1 



2^2 



N 



(4.12) 



This result may be easily generalized to the case of n-fold degeneracy of the extreme slope 
potential curves with an arbitrary n. A simple generalization is possible under conditions 



kj 



CkVl 



j >n , 



k = 1, . . . , 77, , 



(4.13) 



which state that the interaction of degenerate states 1, 2, . . . , n with non-degenerate states 
{j > n) exhibit the same j-pattern, up to common factors Ck- Under these conditions for an 
arbitrary term in (13. ip we obtain 



2p 



^ ^ ttP / ^ 1 \ 



n2 5Z 



Vfe^-^l, 2,...,n "^-^ '^'^ 



-CiiV^uVkxX 1 
p. 



= 



(4.14) 



where = J2k=i ^l- survival probability here we have 
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c- 



N 



exp 



TrVifcVfci 
\/3i-Pk\ 



(4.15) 



fc^l,2,...,n 

We now turn to evaluation of transition probabilities between degenerated states 
1, 2, . . . , n. The expansion terms (I4.14p in fact do not depend on which of degener- 
ate states is initially populated. Formally there is subscript 1 in that indicates initial 
population, but it could be replaced by any j = 2, 3, . . . , n without any other change in 
formulas, except for changing couplings Vifcap-i to V^fc2p-i • 
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However, there is difference in the first term of the perturbative expansion (13. ID that 
exphcitly indicates the initial population. Taking this into account, it is easy to write down 
expression for probabilities of transitions within the submanifold of degenerate states: 

2 



Pi, = ^ 



C4 



exp 



^2 T^VlkVkl \ ^ 



fc^l,2,...,n 



k\ 



j = 2,...,n. (4.16) 



C. Alternative derivation via orthogonalization 

Now we consider alternative scheme of derivation for the case when we have only two 
diabatic potential curves with the same slope Pi = P2 = P and (3 = max^iPk} or P = 
minfc{/3fc}, and ei = 62- As the conditions on couplings we again use formula (14. 6p . 

Introduce a new basis with the states |i) and |2): 

\l) = h{c2\l)-\2)), (4.17) 
\~2) = h{\l) + C2\2)), (4.18) 
h={l + cl)-'^' (4.19) 

instead of states |1) and |2); all other states coincide in the new and old bases. Obviously, 
the new basis is orthonormal. The non-diagonal elements of Hamiltonian matrix with the 
states |1) are all zero: 

{l\H\j)=0, j = 2, 3, 4, ... AT; (4.20) 

in other words state vector |i) is orthogonal to all vectors H\j). This means that the state 
|i) is fully decoupled from all the other states. The diagonal elements of Hamiltonian matrix 
remain the same in new basis. In terms of S matrix this could be written as 

{l\S\i) = 1 , (4.21) 
{2\S\2) = exp (^-n 1(21^^1^)1'^^^^ = V , (4.22) 



where we define 



N ,2 



P^exp -.(l + c^)5:-p^ . (4.23) 
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Here we used the result fl3.23p obtained above for the non-degenerate case. The desired 
S matrix element in the original basis is 

(IIS"!!) = h^cl{l\S\l) + h^{2\S\2) = h\cl + V) . (4.24) 

This gives the state-to-state transition probability 

P,, = h\cl + Vf , (4.25) 

which coincides with the earlier obtained result in (14.121) . 



D. Fully degenerate multistate model 



Consider the case when two fully degenerated bunches of potential curves cross each 
other. The Hamiltonian of this model has form: 

f El ■■■ V ■■■ V\ 



H 



■■■ El V ■■■ V 
V V ■■■ V E2 ■■■ 



(4.26) 



\V V ■■■ V ••• E2/ 

Let n be the number of potential curves with energy Ei = Pit and m be the number of 
potential curves with energy E2 = f32t- The Hamiltonian matrix has dimension {n + m) x 
(n + m). Some transition probabilities for this model can be written down straight off as 
particular cases of formulas (14.15P and fl4.16p . The survival probabilities are 



Pa = At (p"""/^ + m-l)^ j = n + 1, . . . ,n + m . 

The intraband transition probabilities are 

= ^(p"™/'-l)' J^k, j,k = l,...,n, 

Pjk = ^(p'^'^/'-l)' 3^k, j,fc = n + l,...,n + m, 



(4.27) 



(4.28) 
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where p is the standard Landau-Zener probabihty: 

The remaining (interband) probabihties one can obtain by using the normahzation con- 
dition: 

n+m n+m 

E = Y.p^^ = ^- (4.30) 

i=i k=i 

From general considerations it can be concluded that all interband transition probabilities 
are equal, i.e.: 

Pjk = Pjk' , i = 1, ■ ■ ■ , « , k,k' = n + 1, . . . , n + m , 

Pjk = Pjk' , J = n + 1, . . . , n + m , k, k' = 1, . . . , n . (4-31) 

Using fICTD and fOB we obtain: 

P,fe = Pfc, = — , j = l, ...,n, k = n + l, ...,n + m. (4.32) 

nm 

Thus in this highly degenerate multistate model there are only five different state-to-state 
transition probabilities defined by expressions fl4.27p . f l4.28p and (14.321) . This conclusion as 
well as quantitative results were tested by numerical calculations. 



V. CONCLUSION 

In this paper we consider calculation of state-to-state transition probabilities in the gen- 
eralized multistate Landau-Zener model by summation of perturbation theory series. Due 
to specifics of generalized Landau-Zener Hamiltonian (linear growth with time), some of the 
integrals emerging in the pertubative expansions are singular and require special analysis. 
The singularities of these integrals are 'useful' in the sense that they effectively cancel other 
integrations, such that the analytical expressions are obtained for each term in the pertur- 
bative expansion. Subsequently, entire infinite series is summed with the result obtained in 
closed form. The technique of such calculations is one of the principal results of the present 
study. 

The other group of results refers to the degenerate cases. In the general non-degenerate 
case we are able to evaluate only two transition probabilities: the survival probabilities 
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for diabatic potential curves with maximum and minimum slope. Such a situation when 
some state-to-state transition probabilities are expressed by simple analytical formulas, while 
others remain unknown is quite unconventional, although now we know another similar 
example: the multistate Coulomb model [22]. As long as the degeneracy conditions are 
introduced, the analytical expressions for some new state-to-state transition probabilities are 
obtained. For the case of extreme degeneracy, when two fully degenerate bands of diabatic 
potential curves cross each other, the full set of state-to-state transition probabilities was 
derived. Various degenerate cases are met in practice, for example, in the treatment of 



second order effects in Rydberg H atom in perpendicular electric and magnetic fields 21 1. 
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